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I.  INTRODUCTION 


The  Navstar /Global  Positioning  System  is  a satellite -based  navigation 
system  that  provides  extremely  accurate  three-dimensional  position  and 
velocity  information  to  properly  equipped  users  anywhere  on  or  near  the 
earth.  It  is  a Joint  Service  Program,  managed  by  the  Air  Force  with  depu- 
ties from  the  Navy,  Army,  Marines,  Defense  Mapping  Agency,  and  Coast 
Guard  and  with  technical  support  provided  by  The  Aerospace  Corporation. 

The  baseline  program  is  divided  into  three  phases: 

I  - Concept  Validation  Phase  (1974-1979) 

II  - System  Validation  Phase  (1979-1983) 

III  - Production  Phase  (1983-1987) 

The  major  elements  comprising  the  navigation  payload  on  the  satel- 
lites are  the  pseudo  random  noise  signal  assembly  (PRNSA ),  atomic  fre- 
quency standard,  processor,  and  L-band  antenna.  The  PRNSA  includes  the 
baseband  generator,  which  produces  the  basic  P (precise)  and  C/A  (coarse/ 
acquisition)  ranging  codes  and  encodes  navigation  data  from  the  processor 
onto  the  pseudo  random  noise  ranging  signal;  the  amplifier /modulator  units 
that  supply  the  Lj  (1575.42  MHz)  and  Lz  (1227.  6 MHz)  carrier  frequencies 
modulated  by  the  PRN  ranging  signals;  and  the  high-power  amplifiers  that 
amplify  the  carrier  signals  for  transmission. 

A user  measures  pseudo  range  and  pseudo  range  rate  using  the  navi- 
gation signal  from  each  of  four  satellites.  (Pseudo  range  is  the  true  distance 
from  the  satellite  to  the  user  plus  an  offset  due  to  the  user's  clock  bias. 
Similarly,  pseudo  range  rate  is  the  true  slant  range  rate  plus  an  offset  due 
to  the  frequency  of  the  user's  clock.  ) Each  signal  tarries  ephemeris  data 
and  system  timing  information  for  that  satellite,  which  allows  the  user 
receiver /proce ssor  to  convert  the  pseudo  range  and  pseudo  range  rate  to 
user  three-dimensional  position  and  velocity. 
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The  control  segment  consists  of  a Master  Control  Station  (MCS),  a 
navigation  message  upload  station,  and  widely  separated  monitor  stations. 
The  monitor  stations  passively  track  all  satellites  in  view  and  accumulate 
ranging  data,  which  is  processed  at  the  MCS  to  calculate  the  satellite 
ephemerides  and  clock  offsets.  At  least  once  a day  this  information  is 
transmitted  by  the  upload  station  to  the  satellites  for  subsequent  downlink 
transmission  of  the  navigation  data  encoded  on  the  carrier  signals. 

Under  contract  to  the  Space  and  Missile  Systems  Organization  (SAMSO) 
of  the  U.S.  Air  Force,  the  Navstar  satellites  were  developed  and  produced 
by  Rockwell  International  Corp.,  Seal  Beach,  Calif.  Magnavox  Government 
and  Industrial  Electronics  Co.,  Torrance,  Calif.,  has  developed  a variety 
of  user  equipments,  including  the  receivers  now  being  used  at  the  monitor 
stations.  General  Dynamics,  Electronics  Division,  San  Diego,  Calif., 
which  is  responsible  for  the  control  segment,  provided  magnetic  data  tapes 
that  were  used  in  the  preparation  of  this  report. 


II.  SUMMARY 


The  results  of  the  analyses  and  computations  described  in  this  paper 
are  summarized  as  follows: 

a.  The  rms  short-term  error  in  the  pseudo  range  measure- 
ments is  typically  about  0.  3 m.  These  measurements 
are  made  every  6 sec  at  the  monitor  stations,  and  this 
random  error  appears  as  uncorrelated  white  noise.  A 
major  part  of  the  error  is  due  to  the  quantization  in  the 
code  tracking  process. 

b.  When  the  two -frequency  ionospheric  correction  is  added 
to  the  pseudo  range  measurements,  the  error  doubles  to 
about  0.6  m.  This  increase  in  random  error  is  a trivial 
price  to  pay  for  the  overall  improvement  in  accuracy 
associated  with  correcting  for  the  ionospheric  delay. 

c.  The  noise  in  the  delta  range  measurements  over  the  6-sec 
interval  is  about  1 cm  (rms). 

d.  With  a multiplicity  of  satellites  and  monitor  stations,  the 

random  error  in  the  delta  range  measurements  can  be 
separated,  and  the  performance  of  the  atomic  frequency 
standards,  both  in  the  satellite  and  at  the  monitor  stations, 
can  be  evaluated.  At  the  6-sec  sampling  interval,  the 
fractional  frequency  variation  of  these  frequency  standards 
is  about  3 - 4 X which  exceeds  the  specification 

(1  X 1 0 - 1 1 ). 


3 


in.  OVERVIEW  OF  NAVSfAR  PHASE  I ORBIT  CONFIGURATION 


The  baseline  orbit  configuration  for  Navstar  Phase  I is  shown  in 
Table  3-1.  Navstar  satellites  1,  2,  and  3 are  currently  on-orbit.  Navstar  1 
occupies  position  1,  Navstar  2 position  5,  and  Navstar  3 position  6.  For 
convenience,  this  data  is  referenced  to  the  beginning  of  the  day,  midnight 
(0  hr)  Greenwich  Mean  Time  (GMT)  on  1 January  1979.  On  this  particular 
date  the  reference  ground  trace  of  satellite  1 will  first  cross  the  equator 
(north  to  south)  at  a longitude  cf  47. 0 deg  at  4 hr  and  44.  5 min  GMT. 
Approximately  11  hr  and  58  min  later  it  will  again  cross  the  equator  (south 
to  north),  this  time  at  a longitude  of  227.  0 deg.  The  time  of  the  ascending 
node  crossings  will  differ  for  other  days  during  the  year. 

The  ground  trace  of  the  satellite  orbits  is  fixed  from  day  to  day  and 
therefore  repeats  itself  every  23  hr  55  min  and  56.6  sec.  Most  of  the 
deviation  from  a 24-hr  period  is  due  to  the  difference  between  a solar  and 
a siderial  day,  but  a small  part  is  due  to  the  rotation  of  the  orbit  planes 
due  to  the  earth's  oblateness  and  to  sun  and  moon  effects.  This  longitudinal 
motion,  or  precession  effect,  is  slightly  different  for  each  orbit  but  averages 
about  11.5  deg  per  year.  The  satellites  thus  appear  4 min  and  3.4  sec 
earlier  each  day. 

The  information  in  Table  3-1  can  be  used  to  compute  the  times  of 
visibility  of  the  satellites  for  any  location  on  earth  as  of  1 January  1979. 

For  other  days  during  1979,  the  times  are  advanced  4 min  and  3.4  sec  for 
each  day  after  1 January. 


Table  3-1.  Navstar  Phase  I Orbits  at  First  Ascending  Node  on  1 January  1979 


IV.  SIGNAL  STRUCTURE 


Each  satellite  transmits  a navigation  signal  on  two  L-band  frequencies, 
one  at  1575.42  MHz  (L^)  and  the  other  at  1227.6  MHz  (L^).  These  two 
carrier  frequencies  are  biphase  modulated  by  pseudo  random  sequences 
providing  a spread  spectrum  modulation.  The  Lj  carrier  is  actually  modu- 
lated by  two  such  sequences  in  phase  quadrature  so  that,  strictly  speaking, 
this  carrier  is  actually  quadraphase  modulated.  One  pseudo  random  se- 
quence is  a precision  (P)  signal  at  a random  pulse  repetition  rate  of  10.23 
MHz  and  is  an  extremely  long  code  so  that  for  all  practical  purposes  it  is  a 
truly  random  sequence.  The  second  pseudo  random  sequence  is  a coarse 
acquisition  (C/A)  signal,  which  is  a short  sequence  used  either  for  initial 
acquisition  of  the  P signal  or  as  a less  accurate  navigation  signal  for  low- 
cost  users.  The  carrier  frequency  is  biphase  modulated  only  by  the  P 
signal  or,  as  a ground-controlled  option,  only  by  the  C/A  signal. 

For  the  purpose  of  this  discussion  we  will  confine  our  attention  only 
to  the  P signal  on  the  Lj  carrier.  This  is  the  primary  navigation  signal 
used  for  tracking  the  Navstar  satellites  at  the  monitor  stations.  Note  that 
the  carrier  frequency  (1575.42  MHz)  is  an  exact  multiple  (154)  of  the  pseudo 
random  sequence  pulse  rate  (10.23  MHz).  The  wavelength  of  the  carrier  is 
only  19  cm,  whereas  the  "chipping"  rate  of  the  pseudo  random  sequence  is 
about  30  m. 

The  pseudo  random  sequence  is  generated  by  a feedback  shift  register, 
the  output  of  which  modulates  the  carrier  as  illustrated  in  Fig.  4-1.  The 
P signal  pseudo  random  sequence  generator  is  functionally  illustrated  in 
Fig.  4-2.  By  combining  four  12-stage  feedback  shift  registers,  the  equivalen 
of  a 48-stage  shift  register  is  obtained.  Using  a pseudo  random  sequence  to 
biphase  modulate  the  carrier  results  in  the  transmitted  spectrum  being 
spread  as  illustrated  in  Fig.  4-3.  This  spread  permits  an  interference 
signal  from  being  rejected  in  the  receiver  in  the  following  manner:  The 
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receiver's  pseudo  random  sequence  generator  modulates  the  incoming 
signal  in  the  same  manner  as  the  generator  on  the  satellite  (the  two  genera 
tors  are  the  same).  The  original  carrier  frequency  is  thus  reconstri  cted 
and  is  collapsed  to  a very  narrow  band.  The  interference  signal,  however, 
is  spread  out  over  a wide  spectrum,  and  this  signal  can  therefore  be  fil- 
tered out  so  that  only  a small  residue  remains  near  the  now -reconstructed 
carrier  frequency. 
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Fig.  4-1.  Pseudo  Random  (Noise)  Code 
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Fig.  4-3.  Pseudo  Random  Noise-Spread  Spectrum 


V.  GENERIC  RECEIVER 


Functionally,  Navstar  receivers  used  at  the  monitor  stations 
incorporate  two  tracking  loops  that  must  operate  simultaneously  to  properly 
track  the  Navstar  navigational  signal.  The  first  is  the  code  tracking  loop, 
which  tracks  the  pseudo  random  sequence  by  matching  the  locally  generated 
sequence  with  the  sequence  on  the  received  signal.  Simultaneously,  a phase 
lock  loop  is  tracking  the  carrier  frequency.  The  actual  process  is  much 
more  complicated  because  of  the  necessary  intermediate  frequency  (IF) 
downconver sion  steps.  For  simplicity,  however,  these  downconver sion 
steps  are  omitted  in  the  functional  diagram  shown  in  Fig.  5-1.  Figure  5-2 
expands  on  the  function  of  such  a receiver  by  illustrating  the  role  of  a feed- 
back shift  register  and  envelope  detectors  in  the  code  lock  loop  that  tracks 
the  incoming  pseudo  random  sequence.  There  are  three  outputs  of  the  feed- 
back shift  register:  an-on  time  sequence  Pq,  an  early  sequence  P^.,  and 
a late  sequence  Pr  . The  early  and  late  codes  modulate  the  carrier  fre  - 
quency  C,  which  is  synthesized  in  the  phase  lock  loop.  These  signals  are 
then  mixed  with  the  incoming  signal,  thereby  generating  voltages  propor- 
tional to  the  extent  that  the  sequences  match  the  incoming  sequence  from  the 
satellite.  The  difference  between  these  two  signals  generates  an  error 
voltage  that  drives  a voltage  control  oscillator  (VCO)  with  which  the  feedback 
shift  register  is  synchronized,  thereby  tracking  the  incoming  random 
sequence  (Fig.  5-3). 

The  on-time  pseudo  random  sequence  is  mixed  with  the  incoming 
signal  to  reconstruct  the  carrier  signal.  The  phase  lock  loop  (which  is 
actually  a bistable  Costas  loop)  tracks  this  satellite -transmitted  carrier. 
Binary  data  is  added  modulo-2  to  the  P signal  pseudo  random  sequence  at  a 
rate  of  50  bps.  Since  only  the  pseudo  random  sequence  is  removed  from 
the  incoming  signal,  the  50-bps  data  sequence  still  remains  on  the  signal, 
but  the  single -sided  bandwidth  of  this  signal  is  now  only  about  50  Hz. 

In  addition  to  maintaining  phase  lock  on  the  carrier  signal,  the  Costas  loop 
also  strips  off  the  data  remaining  on  the  signal. 
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Fig.  5-1.  Simplified  Diagram  of  Generic  Receiver 
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Fig.  5-2.  Generic  Pseudo  Random  Noise  Receiver 
Functional  Block  Diagram 
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VI.  NAVSTAR  MONITOR  STATIONS 


Each  monitor  station  has  a receiver  that  tracks  both  the  pseudo 
random  sequence  and  the  carrier  signal.  Tracking  the  incoming  pseudo 
random  sequence  allows  the  receiver  to  make  a range  measurement  to  the 
satellite.  If  the  atomic  frequency  standard  at  the  monitor  station  were 
synchronized  exactly  with  the  standard  in  the  satellite,  this  measurement 
would  be  a true  measure  of  the  distance  from  the  monitor  station  to  the 
satellite.  In  a practical  situation,  however,  there  is  always  some  bias 
between  the  monitor  stations  and  the  satellite  atomic  frequency  standards, 
and  because  of  this  relative  clock  difference  the  range  measurement  is 
referred  to  as  pseudo  range.  This  bias  is  determined  as  part  of  the  ephe- 
meris  and  clock  parameter  estimation  process  in  the  MCS  computer,  and 
since  it  is  known  to  some  degree  of  accuracy,  the  pseudo  range  measure- 
ments can  be  considered  a measure  of  the  absolute  distance  from  the  monitor 
station  to  the  satellite. 

Tracking  the  carrier  signal,  however,  does  not  provide  absolute 
ranging  measurements;  only  changes  in  the  range  can  be  measured  from 
carrier  tracking.  At  the  Navstar  monitor  stations,  the  changes  in  range 
are  measured  in  the  phase  lock  loop  over  successive  6-sec  intervals. 

These  range  changes  are  referred  to  as  delta  range  measurements. 

An  important  device  used  in  both  the  code  and  carrier  tracking  loops 
is  an  incremental  phase  modulator  (IPM),  a digital  frequency  synthesizer 
that  is  in  effect  the  actual  VC O.  The  IPM  performs  relative  to  the  frequency 
standard  in  the  monitor  station.  The  quantization  of  each  step  is  1/64  of 
the  wavelength  of  the  quantity  being  tracked.  Since  the  chipping  rate  of  the 
pseudo  random  sequence  corresponds  to  about  30  m,  the  quantization  of  the 
code  tracking  loop  is  about  0.46  m.  Similarly,  since  the  wavelength  of  the 
carrier  frequency  is  about  19  cm,  the  quantization  in  the  phase  lock  loop 
is  about  0.3  cm.  Because  of  this  two -orde r -of -magnitude  difference  in 
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quantization  between  the  code  and  carrier  tracking  loops,  the  tracking 
accuracy  of  the  phase  lock  loop  is  much  greater  than  that  of  the  code  loop. 
The  short-term  quality  of  the  Navstar  tracking  data  is  obtained  by  analysis 
of  the  6-sec  delta  range  data  obtained  from  carrier  tracking  at  the  monitor 
stations . 


VII.  DELTA  RANGE  TRACKING  DATA 
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In  the  laboratory,  atomic  frequency  standards  are  tested  by  comparing 
the  output  of  the  frequency  standard  under  test  with  another  that  (ideally)  is 
much  more  stable.  The  delta  range  tracking  data  makes  possible  the  same 
type  of  test.  The  high-performance  cesium  beam  atomic  frequency  standard 
at  the  monitor  station  serves  as  a reference  for  evaluating  the  atomic  fre- 
quency standard  in  the  satellite  even  though  they  are  about  20,  000  km  apart. 

A quantization  of  0.  3 cm  over  a 6-sec  interval  corresponds  with  an 
average  range  rate  of  0.  5 X 10  m/sec.  Normalized  by  the  velocity  of  light 
(3  X 108  m/sec),  this  corresponds  with  a fractional  frequency  offset  of  the 
atomic  frequency  standard  of  1.7  X 10  ^ “L  Over  6-sec  intervals,  typical 
atomic  frequency  standards  (rubidium  tube  or  cesium  beam)  can  be  expected 
to  exhibit  frequency  stability  on  the  order  of  1 0 ; consequently,  quantization 

of  the  Navstar  tracking  data  is  consistent  with  its  use  for  verifying  proper 
performance  of  the  satellite  atomic  frequency  standard. 


VIII.  POLYNOMIAL  FITTING  PROCESS 


Successive  pseudo  range  and  delta  range  measurements  change  in  a 
gradual  systematic  fashion,  mainly  because  of  satellite  motion  and  earth 
rotation.  There  is  also  some  variation  in  these  measurements  due  to  the 
long-term  effects  of  biases  in  the  satellite  and  monitor  station  frequency 
standard  parameters.  To  evaluate  the  short-term  fluctuations  in  the  fre- 
quency standards,  it  is  necessary  to  remove  these  longer  term  systematic 
effects.  A straightforward  technique  for  accomplishing  this  is  to  fit,  in  a 
least  squares  sense,  a polynomical  function  of  time  to  the  tracking  data. 

The  order  of  the  polynomial  can  vary  from  about  five  to  fifteen,  depending 
on  the  length  of  the  span  of  data  being  used  in  the  analysis. 

A complete,  uninterrupted  span  of  data  can  be  as  long  as  six  hours; 
typically,  however,  good  uninterrupted  spans  are  about  one  to  three  hours. 
The  residuals  of  the  original  tracking  data  to  the  polynomial  fit  are  com- 
puted from  the  polynomial  expressions.  These  residuals  represent  the 
noise  in  the  data  and  can  be  regarded  as  a measure  of  the  short-term  sta- 
bility of  the  frequency  standards.  This  noise  is  generally  similar  to  uncor- 
related white  noise  and  is  not  sensitive  to  the  degree  of  the  polynomial  used 
in  the  fitting  process. 

Since  there  are  two  satellites  and  three  monitor  stations,  by  a covari- 
ance analysis  of  the  residuals  from  various  simultaneous  sets  of  tracking 
data,  it  is  possible  to  separate  the  noise  contributed  by  the  satellite,  that 
contributed  by  the  monitor  station  clocks,  and  other  sources  of  noise 
(usually  the  receiver  tracking  process,  including  the  0.3-cm  quantization). 


IX.  AN  EXAMPLE  OF  TWO -SATE LLITE  TRACKING  DATA 


On  20  July  1978,  there  was  a 70-min  interval  during  which  data  was 
obtained  simultaneously  from  both  Navstars  1 and  2 at  both  the  Vandenberg 
and  Hawaii  monitor  stations.  A complete  covariance  analysis  of  the  resi- 
duals from  a ninth -order  polynomial  fit  to  the  delta  range  data  was  per- 
formed. From  this  analysis,  a breakdown  of  the  residuals  has  been 
obtained.  The  total  noise  from  each  set  of  measurements  was  divided 
between  the  noise  from  the  satellite,  noise  from  the  monitor  stations  clocks, 
and  noise  from  other  sources. 

For  each  of  the  four  sets  of  data,  which  covered  exactly  the  same 
70-min  interval,  7 00  delta  range  measurements  were  analyzed.  The  rms 
of  the  residuals  for  each  set  was  about  1 cm.  The  average  of  the  squares 
of  the  residuals  was  computed,  and  the  average  of  all  six  combinations  of 
the  products  of  the  residuals  was  obtained.  The  computation  can  be  con- 
veniently formulated  in  matrix  form.  The  four  sets  of  700  measurement 
residuals  can  be  regarded  as  a 4 V 700  matrix  A,  an  array  of  4 columns  and 
700  rows.  This  covariance  matrix  is  a 4 X 4 array  C,  which  represents  the 
expected  values  of  the  squares  and  products  of  the  measurement  errors,  where 
the  errors  are  taken  to  be  the  residuals  computed  from  the  ninth-order  poly- 
nomial function  of  time.  The  relationship  between  the  measurement  residuals 
and  the  covariance  matrix  is  given  by 

C = AAt  -r  700 

The  values  of  the  individual  terms  of  the  covariance  matrix  are  computed 

from 


c . . 
U 


21 


1 


For  the  off-diagonal  terms  (the  products  of  different  sets  of  measurements), 
the  normalized  correlation  coefficient  is  given  by 
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For  the  purpose  of  presenting  results,  the  normalized  correlation  co- 
efficients are  sometimes  shown  in  the  upper  right-hand  triangle  of  the  co- 
variance  matrix.  It  should  be  noted  that  the  percentage  correlation  can  be 
either  positive  or  negative  since  the  average  of  the  product  of  two  sets  of 
measurement  errors  can  have  either  sign.  The  quantities  of  interest  are 
frequently  the  one-sigma  of  the  measurement  errors,  and  for  this  reason 
the  square  roots  of  the  variances  along  the  diagonal  are  frequently  presented. 
Similarly,  the  square  roots  of  the  absolutes  of  the  off-diagonal  terms  may  be 
of  interest  and  are  thus  presented  in  the  lower  left-hand  triangle  of  the  modi- 
fied covariance  matrix. 

The  results  of  this  straightforward  covariance  analysis  are  given  in 
Fig.  9-1.  In  this  4x4  modified  covariance  matrix,  the  normalized  correla- 
tion coefficients  in  percent  are  given  in  the  upper  right.  Along  the  diagonal 
of  the  matrix  and  at  the  lower  left,  the  square  roots  of  the  absolute  values  of 
the  elements  of  the  covariance  matrix  are  shown.  The  quantities  on  the 
diagonal  are  the  computed  rms  of  the  residuals  of  the  four  sets  of  measure- 
ments. The  quantities  at  the  lower  left  are  the  parts  of  the  residuals  that  can 
be  attributed  to  isolated  sources.  For  example,  the  0.  553  cm  in  the  first 
column  of  the  second  row  can  be  attributed  to  noise  in  the  cesium  beam  atomic 
frequency  standard  at  Vandenberg.  The  residuals  for  the  cesium  beam  stan- 
dard at  Hawaii  and  the  rubidium  frequency  standards  in  the  two  satellites  were 
also  isolated.  The  cross  correlation  between  the  various  sets  of  data  varies 
from  34.  7 to  44.  1 percent.  Note  that  there  appears  to  be  a 2.  25  percent  cross 
correlation  between  Navstar  1 /Vandenberg  and  Navstar  2/Hawaii  as  well  as 
-0.68  percent  between  Navstar  2/Vandenberg  and  Navstar  1 /Hawaii.  The 
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CROSS  CORRELATION  IN  PERCENT 


0.916 

685 

36.685 

2.251 

NAVSTAR  1/ VANDENBERG 

0.553 

o7%4^> 

^-0.683 

40.244 

NAVSTAR  21  VANDENBERG 

0.591 

0.083 

1.038^ 

■^44.089 

NAVSTAR  1/ HAWAII 

0.152 

0.658 

0.715 

1.  H7  " 

NAVSTAR  2IHAWAII 

NAVSTAR  1/ 
VANDENBERG 

NAVSTAR  21 
VANDENBERG 

NAVSTAR  11 
HAWAII 

NAVSTAR  21 
HAWAII 

Fig.  9-1.  Modified  Covariance  Matrix 

expected  cross  correlation  for  these  cases  is,  of  course,  zero,  and  the  fact 
that  finite  values  were  computed  simply  reflects  the  statistical  variability  of 
a limited  set  of  observation  data. 

The  results  shown  in  Fig.  9-1  are  presented  in  a different  form  in 
Table  9-1.  Noise  from  various  atomic  frequency  standards  and  other  sources 
are  separable  when  there  is  a multiplicity  of  satellites  and  monitor  stations. 
These  results  show  that  at  a 6-sec  sampling  interval,  all  four  frequency  stan- 
dards have  comparable  stability,  about  3-4X10  ^ fractional  frequency 

variation.  The  category  "other"  was  obtained  by  subtracting,  in  an  rms 
sense,  the  satellite  and  monitor  station  clock  residuals  from  the  total  residual. 

The  observed  fractional  frequency  stability  over  a 6- sec  sample  time 
exceeds  the  specification  level  of  10  ^ and  is  typical  of  the  performance 
measured  by  Rockwell  International  in  the  laboratory.  In  Fig.  9-2,  based 
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Table  9-1.  Breakdown  of  Residuals  (rms, 
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Fig.  9-2.  Frequency  Standard  Stability  Performance 

on  Rockwell  data,  the  specification  for  the  satellite  rubidium  frequency 
standard  is  shown  as  a heavy  solid  line.  The  many  dots  represent  the 
results  of  laboratory  test  data  taken  from  seven  frequency  standards.  Note 
that  at  or  near  a 6 -sec  sample  time,  the  data  points  cluster  in  the  region 
from  3 - 6 X 10’  . In  conclusion,  the  analysis  of  the  satellite  tracking 

data  indicates  substantially  the  same  level  of  frequency  standard  perfor- 
mance as  generally  observed  in  the  laboratory. 

Figure  9-3  shows  a sample  of  the  residuals  for  a 10-min  interval 
over  which  100  measurements  were  taken.  This  plot  represents  the  first 
10  min  of  the  data  taken  at  the  Vandenberg  monitor  station  tracking  Navstar  1. 
This  plot  is  typical  of  all  the  data  obtained  in  tracking  both  satellites  from 
both  monitor  stations. 
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RANGE  in  ni.IERS 


X.  PSEUDO-RANGE  TRACKING  DATA 


While  the  monitor  station  is  tracking  both  pseudo  range  and  delta 
range  on  the  frequency  (1575.42  MHz),  it  is  also  measuring  the  difference 
between  the  pseudo  range  as  measured  on  the  L ^ and  L ^ (1227. 6 MHz) 
frequencies.  This  difference  is  measured  to  automatically  correct  for  the 
ionospheric  delay  in  the  transmission  of  the  navigation  signal.  The  primary 
reason  that  two  frequencies  are  employed  on  Navstar  is  to  make  possible 
this  compensation,  which  is  based  on  the  fact  that  the  ionospheric  delay  is 
inversely  proportional  to  the  square  of  the  carrier  frequency.  The  dif- 
ference in  the  delay  as  measured  on  the  two  L-band  frequencies  is,  in  effect, 
a measurement  of  the  ionospheric  delay  at  each  of  these  two  frequencies. 

The  expression  for  the  necessary  compensation  of  the  pseudo  range  mea- 
surement at  the  Lj  frequency  is  given  by 


or 

R'  = R + 1.  54573  (R(  - R2) 

whe  re 

R ' - corrected  pseudo  range  measurement  at  the 
Lj  frequency 

Rj,  R£  = receiver  outputs,  i.e.,  uncorrected  pseudo  range 
measurements  at  the  Lj,  frequencies 

The  short-term  quality  of  the  pseudo  range  measurements  can  be 
obtained  by  the  same  technique  previously  described  for  the  delta  range 
measurements.  The  quality  can  be  evaluated  for  both  the  uncompensated 
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pseudo  range  measurement  and  the  pseudo  range  measurement  corrected 
for  the  ionospere  R'.  From  this  expression,  the  noise  in  the  compensated 
pseudo  range  measurements  would  be  expected  to  be  larger  than  the  mea- 
surements before  the  ionospheric  compensation  because  there  are  two  noise 
sources,  the  noise  in  tracking  the  pseudo  random  code  on  the  L{  frequency 
and  the  noise  in  the  receiver  tracking  of  the  difference  in  pseudo  range  on 
the  two  L -band  carrier  frequencies. 

For  the  7 0-min  interval  on  20  July  1978,  a ninth -order  polynomial 
function  of  time  was  fitted  to  the  pseudo  range  data  taken  at  the  Vandenberg 
monitor  station  while  tracking  Navstar  1.  For  the  uncompensated  pseudo 
range  measurements  the  rms  of  the  residuals  was  computed  to  be  0.29  m, 
and  for  the  ionospheric -compensated  pseudo  range  measurements  the  rms 
was  computed  to  be  0,57  m.  These  results  represent  both  the  code  tracking 
of  the  pseudo  random  sequences  and  the  noise  increase  associated  with 
adding  the  ionospheric  correction.  A representative  10-min  sample  of  these 
residuals,  which  covers  1 00  measurements,  is  shown  in  Figs,  10-1  and 
10-2  for  the  uncompensated  and  compensated  cases,  respectively. 
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Fig.  10-1.  Uncompensated  Pseudo  Range  Residuals 
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Fig.  10-2.  Compensated  Pseudo  Range  Residuals 


PSUEDO  RANGE  IN  METERS 


XI.  ANOTHER  EXAMPLE  OF  TWO 
SATELLITE  TRACKING  DATA 


On  22  August  1978,  there  was  a two -hour  span  during  which  data  was 
obtained  simultaneously  from  both  Navstars  1 and  2 at  all  three  monitor 
stations,  with  high-quality  data  for  all  six  combinations.  The  residuals 
from  a ninth -order  polynomial  fit  of  the  delta  range  data  (la)  are  given  in 
Table  11-1. 

Table  11-1.  Residuals  from  a Ninth-Order  Fit 


Combination 


Navstar  1/Vandenberg 
Navstar  2/Vandenberg 
Navstar  1 /Hawaii 
Navstar  2/Hawaii 
Navstar  1 /Alaska 
Navstar  2/Alaska 


Residuals  (la). 


0.  976 


0.  960 
1. 039 
1. 097 
0.  950 


1. 055 


A covariance  analysis  of  these  residuals  was  performed,  and  the  results 
are  given  in  Fig.  11-1.  As  in  the  previous  example,  the  6x6  matrix 
shown  on  this  figure  is  a modified  covariance  matrix  in  which  the  normalized 
correlation  coefficients  in  percent  are  given  in  the  upper  right-hand  trian- 
gular part  of  the  matrix.  Along  the  diagonal  of  the  matrix  and  the  lower 
left-hand  triangle,  the  square  roots  of  the  elements  of  the  covariance  matrix 
are  shown.  The  percentage  cross  correlations,  the  one -sigma  measurement 
residuals,  and  the  manner  in  which  the  noise  error  can  be  assigned  to  satel- 
lites, monitor  station  clocks,  and  other  sources  is  similar  to  the  results 
obtained  from  the  first  example. 
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MS  OF  RESIDUALS 
IN  CENTIMETERS 


CROSS  CORRELATION  IN  PERCENT 
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Fig.  11-1.  Modified  Covariance  Matrix 


